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We investigate the critical behavior of three-dimensional antiferromagnetic CP^ -1 (ACP^^ 1 ) 
models in cubic lattices, which are characterized by a global U (N) symmetry and a local U(l) 
gauge symmetry. Assuming that critical fluctuations are associated with a staggered gauge-invariant 
(hermitian traceless matrix) order parameter, we determine the corresponding Landau-Ginzburg- 
Wilson (LGW) model. For N = 3 this mapping allows us to conclude that the three-component 
ACP 2 model undergoes a continuous transition that belongs to the 0(8) vector universality class, 
with an effective enlargement of the symmetry at the critical point. This prediction is confirmed by 
a detailed numerical comparison of finite-size data for the ACP 2 and the 0(8) vector models. We 
also present a renormalization-group (RG) analysis of the LGW theories for N > 4. We compute 
perturbative series in two different renormalization schemes and analyze the corresponding RG flow. 

We do not find stable fixed points that can be associated with continuous transitions. 

PACS numbers: 05.70.Jk,05.30.Cc 


I. INTRODUCTION 

Qpjv-i moc i e i s are a c l ass of models in which the 
fundamental field is a complex ./V-component unit vec¬ 
tor (more precisely, an element of the complex projec¬ 
tive manifold CP^ -1 ), and which are characterized by 
a global V(N) symmetry and a local U(l) gauge sym¬ 
metry. They emerge as effective theories of SU(iV) quan¬ 
tum antiferromagnets and of scalar electrodynamics 

with a compact U(l) gauge group. The simplest three- 
dimensional (3D) CP^ -1 lattice model is defined by the 
Hamiltonian 

H = jJ2\z i -z j \ 2 , ( 1 ) 

< ij) 

where the sum is over the nearest-neighbor sites of a cu¬ 
bic lattice, Zi are ./V-component complex vectors satisfy¬ 
ing Zi ■ Zi = 1. The model is ferromagnetic for J < 0, 
antiferromagnetic for J > 0. 

The CP 1 model can be mapped onto the 0(3)- 
symmetric Heisenberg model. Indeed, if one defines 0(3) 
spins sf = Y^ab *i a ab z i ) where a a are the Pauli matrices, 
one can rewrite the CP 1 Hamiltonian as that of the usual 
3-vector Heisenberg model. As a consequence, the critical 
properties can be straightforwardly derived by using the 
wealth of results available for the Heisenberg model, see, 
e.g., Refs. (HQ. On the other hand, several aspects of 
the phase behavior of CP^” 1 models with N > 2 remain 
unclear and worth being further investigated. 

The critical behavior of these models can be inves¬ 
tigated by constructing an effective Landau-Ginzburg- 
Wilson (LGW) theory. This approach requires the iden¬ 
tification of the order parameter associated with the crit¬ 
ical modes. Assuming that the critical behavior is essen¬ 
tially driven by gauge-invariant modes, a plausible choice 
for ferromagnetic (J < 0) systems is the gauge-invariant 


site variable am 

Qt = nA - ( 2 ) 

which is a hermitian and traceless N x N matrix. In 
the corresponding LGW theory, the fundamental field 
is therefore the most general traceless hermitian matrix 
$ ab (a:), which one can imagine being defined as the av¬ 
erage of Q“ b over a large but finite lattice domain. The 
Hamiltonian is the most general fourth-order polynomial 
in $ consistent with the U(-ZV) symmetry: 

H = Tr^T) 2 + r Tr $ 2 + w tr <h 3 (3) 

+ u (Tr $ 2 ) 2 + v Tr $ 4 . 

For N = 2 one recovers the 0(3)-symmetric LGW the¬ 
ory because the cubic term vanishes and the two quartic 
terms are equivalent [§J, consistently with the equivalence 
between the CP 1 and the Heisenberg model. Because of 
the presence of the cubic term, on the basis of mean- 
field arguments, one expects the system to undergo a 
first-order transition for any N > 2, unless the Hamilto¬ 
nian parameters are tuned so that w = 0 in the effective 
model. This prediction is, however, contradicted by re¬ 
cent numerical studies msniEi , which find evidence 
of continuous transitions in models that are expected 
to be in the same universality class as that of the 3D 
CP 2 model. In particular, a numerical study of 3D loop 
models P provided the estimate v = 0.536(13) for the 
correlation-length critical exponent. These results imply 
the existence of a 3D CP 2 universality class characterized 
by a U(3) global symmetry and U(l) gauge invariance, 
with a corresponding fixed point (FP) that cannot be de¬ 
termined in perturbation theory at fixed N. In order to 
access this FP, the authors of Ref. @] proposed a double 
expansion around N = 2 (where the cubic term vanishes) 
and e = 4 — d, arguing that a continuous transition may 
be possible for values of N sufficiently close to N = 2. 
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For larger values of TV, i.e. TV > 4, the numerical analy¬ 
ses i I 0 show evidence of first-order transitions. In 
the large-TV limit, the quantum field theory correspond¬ 
ing to the ferromagnetic CP ^ 1 model gives rise to an 
effective abelian Higgs model and Landau-Ginzburg (LG) 
theory of superconductivity 0 , whose renormalization- 
group (RG) flow presents a stable FP for a sufficiently 
large number of components. Thus continuous transi¬ 
tions are again possible in CP ^ -1 models at large TV, 
sharing the same universal critical behaviors of the LG 
theory of superconductivity 0. These results are again 
in contrast with the conclusions obtained from the LGW 
theory ([3]) and suggest that, at least for large values of 
TV, critical modes are not exclusively associated with the 
gauge-invariant order parameter Q, cf. Eq. 0 . 

In this paper we investigate the critical behavior of an¬ 
tiferromagnetic CP W-1 (ACP^ -1 ) models, such as those 
described by the Hamiltonian 0 with J > 0, on a cu¬ 
bic lattice (we expect a similar behavior on any bipartite 
lattice). For TV = 2, they undergo a critical transition 
in the same universality class as that of the ferromag¬ 
netic CP 1 model. Indeed, the ACP 1 model is equiva¬ 
lent to the antiferromagnetic Heisenberg model, which in 
turn can be mapped onto the ferromagnetic one by per¬ 
forming the transformation s(x) (—l) a:i+:C 2 +a: 3 s(a;), 
where x = (aq, £ 2 , £ 3 ). Therefore, the staggered variables 

Sstag = (-1 )*1+Z2+*3 S or Q stag = (_l)Xl+X2+ 3: 3 Q have 

the same critical behavior as s or Q in the ferromagnetic 
model. However, for TV > 2 the behavior of ACP ^” 1 
models differs from that of ferromagnetic CP ^ -1 mod¬ 
els, as we shall show. 

Keeping the assumption that the critical modes can be 
represented by local staggered gauge-invariant variables, 
we show that the LGW Hamiltonian describing the be¬ 
havior of the critical modes in the ACP ^ -1 models is the 
one given in Eq. 0 , without the cubic term, that is with 
w = 0. Indeed, the staggered nature of the order parame¬ 
ter gives rise to a symmetry <f> —> — <f>, which prevents the 
presence of odd terms in <F, such as the cubic term. This 
fact greatly simplifies the RG analysis of the theory. In 
particular, it allows us to predict that the critical behav¬ 
ior of the ACP 2 model belongs to the universality class 
of the 0 ( 8 ) vector model, with a dynamical enlargement 
of the symmetry at the critical point. Correspondingly, 
we predict v « 0.85 and r\ « 0.03 for the ACP 2 model, 
which differ from those found in Refs. i0 for the fer¬ 
romagnetic CP 2 universality class. To validate this pre¬ 
diction, we perform Monte Carlo (MC) simulations of the 
lattice ACP 2 and 0(8) vector models. The critical expo¬ 
nents and finite-size scaling (FSS) functions turn out to 
be the same in the two models, in agreement with the 
RG argument. Finally, we present a general RG study 
of the LGW theory 0 without cubic term. We com¬ 
pute high-order field-theoretical (FT) perturbative series 
in two different schemes. The RG analysis does not pro¬ 
vide evidence of the existence of stable FP for TV > 4. 

We mention that a critical-point symmetry enlarge¬ 
ment analogous to that of ACP 2 models occurs in the 


antiferromagnetic RP 2 model 00 - These systems 
are similar to those considered here: their Hamiltonian 
is also given by Eq. 0, but the site variable is a real 
unit vector. In the RP 2 case, an analogous RG argument 
based on the corresponding LGW theory shows that the 
model should have a critical behavior in the universality 
class of the 0(5) vector model. 

The paper is organized as follows. In Sec. [TT] we con¬ 
struct the LGW theory which is expected to describe 
the critical modes at continuous transitions of ACP ^ -1 
models, assuming a staggered gauge-invariant order pa¬ 
rameter. Sec. m is devoted to a numerical study of the 
AGP 2 model. We show that its continuous transition is 
in the same universality class as that of the 0 ( 8 ) vec¬ 
tor model. In Sec.|lV]we study the RG flow relevant for 
models with more components, i.e., for TV > 4, comput¬ 
ing and analyzing high-order FT perturbative series for 
the corresponding LGW theories. Finally in Sec. |V] we 
summarize our main results and draw some conclusions. 
Some details are reported in the appendices. 

II. LGW THEORY FOR THE ACP ^" 1 MODELS 

We now derive the LGW theory for the critical modes 
at the transition in ACP ^ -1 models, which is completely 
specified by the symmetry of the model, the nature of the 
order parameter, and the symmetry-breaking pattern. In 
the case at hand, the model has a global U(TV) symmetry 
and a local U(l) gauge invariance. We assume that the 
critical modes are effectively represented by local gauge- 
invariant variables, such as 0. In the case of antifer¬ 
romagnetic interactions (J > 0), the minimum of the 
Hamiltonian 0 is locally realized by taking Zi • Zj = 0 
for any pair of nearest-neighbor sites. 

In order to construct the LGW Hamiltonian, we should 
identify the order parameter of the transition. At vari¬ 
ance with the ferromagnetic case, in the ACP ^ -1 model 
we should take into account the explicit breaking of trans¬ 
lational invariance in the low-temperature phase. To 
clarify the issue, let us consider the antiferromagnetic 
O(Af) vector model with Hamiltonian Ho = Xqij) s i'Sj- 
In this case the order parameter is cf> = ^2, x p x s x , w h ere 
p x is the parity of the site x = (aq, £ 2 , £ 3 ) defined by 
p x = (— l)2Tk Xi . Under translations of one site, we find 
4> —> —</>, hence this parameter allows us to probe the 
breaking of translational invariance. In the ACP ^ -1 
model the natural field variable is the combination zfz b , 
which is invariant under the local U(l) gauge transforma¬ 
tions of the model. Hence, we define the order parameter 

B ab = Y,P,z a x z b x . (4) 

(x) 

It is immediate to verify that B is hermitian and trace¬ 
less [this follows from the presence of p x ], that it changes 
sign under translations of one site which exchange the two 
sublattices, and that it coincides with the 0(3) order pa¬ 
rameter for N = 2. Then, as usual, in order to construct 
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the LGW model, we replace B with a local variable $ as 
fundamental variable (essentially, one may imagine that 
W is defined as B , but now the summation extends only 
over a large, but finite, cubic sublattice). Then, we write 
down the most general fourth-order polynomial that is 
invariant under U (N) transformations and under the Z 2 
transformation T —> — T, a consequence of the transla¬ 
tion invariance of the original theory. We obtain 

H a =Tr(a A1 ^) 2 +rTr^ 2 + ^(TrT 2 ) 2 + ^TrT 4 . (5) 

The original U (N) symmetry corresponds to the symme¬ 
try T —> U'^U where U £ U (N). The order parameter 

is a hermitean traceless matrix as the variable $ intro¬ 
duced in the ferromagnetic case. However, because of the 
presence of the symmetry —> — T, which is a specific 
feature of the antiferromagnetic model, the LGW Hamil¬ 
tonian ([5J does not present a cubic term, which instead 
appears in the <f> 4 Hamiltonian 0 corresponding to the 
ferromagnetic case. 

Note that the model is not only characterized by the 
symmetry group, but also by the nature of the order 
parameter. There are indeed other models with U(iV) 
symmetry (see, e.g., Ref. [h| for an example), which, 
however, have a different order parameter and different 
symmetry-breaking patterns, leading to different univer¬ 
sality classes. 

The stability domain of H a can be determined by 
studying the asymptotic behavior of the potential 

R(«I») = r Tr 'F 2 + ^ (Tr T 2 ) 2 + ^ Tr ^ 4 . (6) 

This analysis can be easily performed by noting that 
E(’F) only depends on the N real eigenvalues A 0 of the 
hermitian matrix T, which satisfy the condition E a A a = 
0. We find that the theory is stable if 

, n , N 2 — 3N + 3 
Uo + OnVo > 0 , b N = jyny- _ ^ , (7) 

and if 

Mo + jf y o > 0 for even N, (8) 

uo + cnv 0 > 0 for odd N, 

where 

N 2 + 3 

CN N{N 2 ~ 1)' ^ 

Physical systems corresponding to the effective theory 
(ED with Uo , Uo that do not satisfy these constraints are 
expected to undergo a first-order phase transition. 

The analysis of the minima of the potential P(1I/) for 
r < 0 gives us information on the symmetry-breaking 
patterns. For vq < 0, the absolute minimum of P(ff') is 
realized by configurations with T = t/T m i n W and 


where I n indicates the n x n identity matrix. This gives 
rise to the symmetry-breaking pattern 

U(1V) U(l) x U(7V - 1). (11) 

On the other hand, for vq > 0 and even N the minumum 
is realized by 


^ min ~ ( -I N , 2 ) (12) 

implying the symmetry-breaking pattern 

U(1V) ^ U(iV/2) x U(lV/2). (13) 

For uo > 0 and odd values of N , we have instead 


^min 


I(N+l)/2 0 

0 —kI(N- 1)/2 


k = (N + 1)/(N — 1), so that 


U(iV) ->• U(iV/2 + 1/2) x U(lV/2 - 1/2). 


(14) 


(15) 


Note that for N = 3, the symmetry breaking patterns 
HID and m are equivalent, hence the sign of vo does 
not play any role. 

An important remark is in order. The derivation of 
the LGW Hamiltonian 0 is based on the assumption 
that the order parameter is the staggered and traceless 
hermitian matrix 0. This assumption can be checked 
for N = 3. As we show in App. E the minimum-energy 
configurations of Hamiltonian 0 for J > 0 have a very 
simple structure. Modulo a global U(3) transformation, 
one can take Zj = ( 1 , 0 , 0 ) on one sublattice, and z,; = 
(0, 0 ^, 6 ,;) on the other one. Hence, a zero-temperature 
configuration corresponds to 


10 0 \ 

b = 0 -E.M 2 -EXfc ■ (is) 

Vo -E i<Hb* -EiN 2 / 


Therefore, B is nonvanishing in the low-temperature 
phase and represents the correct order parameter. The 
symmetry-breaking pattern is therefore that given in 
Eq. (fill) or, equivalently, Eq. (fl5l) . 

For N > 4, we have not been able to identify ordered 
zero-temperature configurations that are translation in¬ 
variant at least on one sublattice, hence, we have not been 
able to check that B, as defined in Eq. 0 , is nonvanish¬ 
ing in the low-temperature phase, hence that it can be 
taken as the order parameter. In the following, we make 
the working hypothesis that this is the case, determin¬ 
ing what this assumption implies for the nature of the 
transitions in the ACP ^ -1 models. 

For N = 3 the LGW theory 0 simplifies. Indeed, one 
can easily prove that 


^min 


In- 1 0 A 

0 -(IV-1) )’ 


( 10 ) 


TrT 4 = i(TrT 2 ) 2 


(17) 
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for any 3x3 traceless hermitean matrix. Then, let us 
define an 8 -component real vector field as follows: 


Tn + a+fI/22 

7f 

Tl2 = <t>3 + i<t> 4 


Til + «-T22 


— <Pl> 


72 

T 13 = fe + iq. 


= <h, ( 18 ) 

T 23 = </>7 + 


where a± = (1 ± 73)/2. In terms of the field <j> we have 


-TrT 2 = 


The, we can rewrite H a as 


Ho = (77) 2 


+ 2r (ft 2 + go (4> 2 ) 2 , 


(19) 


( 20 ) 


where go = uo + vo/2, proving that the model is equiva¬ 
lent to the 0 ( 8 ) vector theory. 

Since the 0(8) vector theory has a stable Wilson-Fisher 
FP, the above correspondence allows us to predict that 
the critical behavior of the ACP 2 model must share the 
same universal features as the 3D 0(8) vector model. 
Note that the enlargement of the symmetry to 0(8) is 
a feature of the LGW theory (0, i.e., of the expansion 
up to fourth powers of T. Indeed, one can easily check 
that the sixth-order terms allowed by the U(3) symmetry, 
such as TrT 6 , do not share the 0(8) symmetry. Since 
these terms are RG irrelevant at the 0(8) FP, the con¬ 
tribution of the terms breaking the 0 ( 8 ) symmetry is 
suppressed at the critical point. In this sense we have 
a dynamic enlargement of the symmetry to 0 ( 8 ) at the 
critical point. 

Of course, the fact that only one quartic term is in¬ 
dependent for TV = 3 holds also for the Hamiltonian ([3]) 
corresponding to the ferromagnetic CP 2 model. Thus, it 
has an 0(8) FP, which is unstable due to the presence 
of the cubic term. This cubic term gives rise to a spin-3 
perturbation of the 0(8) model 17]. Its RG dimension 
2/3 can be estimated using the results of Ref. [l 6 j for the 
crossover exponent </> 3 , i.e. 2/3 = <j> 3 /u. The crossover ex¬ 
ponent </>3 at the 0(8) FP was estimated in Ref. [l(|, ob¬ 
taining <j >3 = 0.97(3) in the fixed-dimension massive-zero- 
momentum scheme, and cf >3 = 0.95(5) in the e expansion. 
Since v = 0.85(2) (see Sec. lIIIBl) . we obtain y 3 « 1.1. As 
?/3 > 0, the 0(8) FP is unstable in the presence of the 
cubic term, and, therefore, the 0(8) FP cannot be the 
stable FP in the case of ferromagnetic interactions. This 
is confirmed by the results of Refs. Sill. The critical 
exponents, for instance v = 0.536(13), significantly differ 
from those of the 0 ( 8 ) universality class. 

It is worth mentioning that an analogous enlarge¬ 
ment of the symmetry at the critical point is also ob¬ 
served 01 " in the antiferromagnetic RP 2 model. The 
Hamiltonian of RP ^^ 1 models is analogous to that of 
CpN-i models It is given by Eq. ©, with real TV- 
component spins Si replacing the complex vectors Zi. 
The order parameter should be a symmetric and traceless 
TV x TV matrix E, analogous to the matrix B defined in 
Eq. ©. The corresponding LGW Hamiltonian is given 
in Eq. ©, with T replaced by E [rJ [b|. For TV = 3, 


the two quartic terms are proportional, and one obtains 
the <f > 4 theory of the 0(5) vector model (again high-order 
terms break the 0(5) symmetry, but since they are irrel¬ 
evant, the 0(5) symmetry breaking is suppressed in the 
critical limit). The numerical results of Refs. [l3|, |l4| for 
the antiferromagnetic RP 2 model confirm this prediction. 

To identify the nature of a critical transition for TV > 4, 
if it exists, we must determine the fixed points of the RG 
flow of the LGW theory ©. The absence of a stable 
FP implies the absence of continuous transitions. Such 
an analysis is quite complex, due to the presence of two 
quartic terms. We will perform it in Sec. m 


III. CRITICAL BEHAVIOR OF THE ACP 2 
LATTICE MODEL 


We now check the predictions of the previous section 
for the critical behavior of the 3D ACP 2 model, confirm¬ 
ing that it undergoes a continuous transition in the 0 ( 8 ) 
universality class. 


A. Monte Carlo simulations and observables 


In order to study the critical behavior of the ACP 2 
lattice model © with J = 1, we perform Monte Carlo 
simulations of cubic systems of linear size L with periodic 
boundary conditions. Because of the antiferromagnetic 
nature of the model we take L even. We use a standard 
Metropolis algorithm, hence we are only able to obtain 
reliable results up to L = 40. We use a simple updating 
algorithm. If ip = (Rez,Im.z) is a six-component vec¬ 
tor, the update consists in proposing the new vector R<p, 
where R is a random 0(2) matrix acting on two randomly 
chosen components of <p. 

In our MG simulations we compute correlations of the 
gauge invariant operator 

K b = ( 21 ) 


Its two-point correlation function is defined as 


G(x — y) = (TrPjP y ) = <|i 


7 


( 22 ) 


Due to the staggered nature of the ordered parameter, we 
should distinguish correlations between points belonging 
to the same sublattice and points belonging to different 
sublattices. Here we define the susceptibility and the 
correlation length by summing only over points with the 
same parity: 


X 


= £ G{x) = G{ 0), 


(23) 


7 = 


G{ 0) - G(p) 


4sin 2 (p min /2) G(p) 


where x runs over all even points, 


G(P)= E 


0 zp-x 


G(x) 


(24) 


(25) 
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FIG. 1: (Color online) MC data of R £ for the ACP 2 lattice 
model for different lattice sizes L. They show a crossing point 
at p « 4.14. The dotted lines are drawn to guide the eye. The 
vertical dashed line corresponds to our best estimate m ° f 
Pc obtained by a FSS analysis of the available data assuming 
the transition to belong to the 0(8) universality class. 

is the Fourier transform of G{x ) over the even sublattice, 
p = (pmin,0, 0), and Pmin = 27 t/L. Finally, we consider 
the Binder parameter 

T T _ (Ex even TrPjPx] 2 } 

cr I r / * /' ' 1 

\2-^x even ^0 x i 

In the FSS limit any RG invariant quantity R, such as 
R^ = £/L and U, is expected to behave as 

R(J3,L) = f R (X) + L~ u g R (X) + ... (27) 

where X = ((3 — /3 C )L 1 ' v and f r{X ) is a universal func¬ 
tion apart from a trivial normalization of the argument. 
In particular, the quantity R* = fn(0) is universal within 
the given universality class. The approach to the asymp¬ 
totic behavior is controlled by the universal exponent 
u> > 0 , which is associated with the leading irrelevant 
RG operator. Around /3 C one may expand fn{X) and 
gn(X) in powers of the scaling variable X, obtaining 

R = R*+J2 bn(p-Pc) n L nl "+L-“ £ c n (p-p c ) n L n ' v +... 

71=1 71=0 

(28) 

The exponent 77 is determined by analyzing the FSS be¬ 
havior of the susceptibility 


FIG. 2: (Color online) FSS behavior of R £ for the 0(8) vector 
model. Plot of R £ versus (/3 — f3 c )L 1 R, setting fi c = 1.92677 
and v = 0.85. The dashed horizontal line corresponds to 
Rs = 0.5237. 

v that gives the slope of the data at the critical point. At 
a first-order transition we expect v = 1/d = l/2> prUrlj . 
while v > 1/d at continuous transitions. The numeri¬ 
cal data, including those of the specific heat, definitely 
exclude a first-order transition. Indeed, the increase of 
dR/d/3 at the crossing point is much slower than L 3 . 
Therefore, we conclude that the ACP 2 model has a con¬ 
tinuous transition. 


B. The 3D 0(8) vector model 

To show that the ACP 2 lattice model belongs to the 
0(8) universality class, we show that critical exponents 
and FSS curves are the same in the two models. Thus, 
we begin by computing these quantities in the 0(8) spin 
model defined by the Hamiltonian 

Hp = ~Y S >- Sj, (30) 

<b) 

where the spin variable Si is an 8-component unit vector. 
We consider cubic systems of linear size L with periodic 
boundary conditions. 

We consider the two-point function 

G 0 {% 7/) (&x ' (31) 


X ~ L 2 ~ v {f x {X) + 0(L-“)] . (29) 

We present a FSS analysis of the numerical data of the 
ACP 2 lattice model, up to L = 40. In Fig. |I] we show 
= £/L for several values of L. Data have a crossing 
point, providing evidence of a transition in the interval 
4-1 P < 4.2. An analogous behavior is shown by the 
Binder parameter U. 

To determine whether the transition is continuous or 
of first order, we should estimate the effective exponent 


and compute the corresponding susceptibility and 
second-moment correlation length. They are defined as 
in Eqs. (E3l) and EH) , but now we sum over all lattice 
points, as the model is ferromagnetic. Moreover, we con¬ 
sider the Binder parameter U defined as in Eq. (ESD : we 
replace P x with s x and sum over all lattice points. 

We perform simulations on lattices of size up to L = 96 
(we use a cluster algorithm) and estimate R^, U, and 
X- To compute the critical exponents we fit the data to 
Eqs. E71) and EH) . In the analysis we take into account 
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the scaling corrections of order L _w , fixing ui ~ 0.8, as 
predicted by the FT perturbative analyses discussed be¬ 
low. Corrections turn out to be small, hence the analyses 
of the MC data are unable to provide a more accurate 
estimate of w. Fits of R^ to Eq. (l2Sl) (we take the first 
terms in the expansions) give /3 C = 1.92677(2), 

^ = 0.85(2), i] = 0.0276(5), (32) 

and the universal critical value i?| = 0.5237(4). The 
quoted uncertainty includes the statistical error and the 
variation of the estimates as u> varies between 0.75 and 
0.85, an interval that is larger than that obtained in the 
FT analyses reported below. In Fig. [5] we show R ^ versus 
(/3 — p c )L using the above-reported estimates. Data 
collapse onto a single curve, confirming the accuracy of 
the estimates. Consistent results are obtained from the 
analysis of the Binder parameter, which also gives U* = 
1.0383(3). 

The values of the critical exponents can be compared 
with previous results. Field theory gives v ss 0.830 and 
r\ ss 0.027 [23], while the analysis of strong-coupling ex¬ 
pansions gives v ss 0.84, 0.86 [23j (the two estimates are 
obtained by means of two different resummation meth¬ 
ods). Within errors, they agree with the estimates (1321) . 
We have repeated the analysis of the available six-loop se¬ 
ries within the massive zero-momentum renormalization 
scheme [13, H3|, using the conformal mapping method 
that exploits the known large-order behavior of the per¬ 
turbative expansions [25!, 26). We obtain v = 0.826(4), 
rj = 0.025(1), and w = 0.81(1), where the errors are re¬ 
lated to the change of the estimates with respect to a 
(reasonable) variation of the parameters entering the re¬ 
summation procedure. They are substantially consistent 
with our favorite MC estimates (l32l) . although one may 
suspect that errors are slightly underestimated. The FT 
analysis also provides the estimate of u> that we used 
(note that, to be on the safe side, we allowed for a much 
larger uncertainty in the MC analysis). 

C. FSS of the ACP 2 lattice model 

If the transitions in the ACP 2 and 0(8) models belong 
to the same universality class, critical exponents and FSS 
curves fn(X) for RG invariant quantities (apart from a 
trivial rescaling of the variable X ) should be the same. 
This is what we check below. 

A first unbiased universality check which does not need 
an estimate of the critical point /3 C , is obtained by plot¬ 
ting versus U. Indeed, since both R^ and U satisfy 
Eq. (EH), we must have 

R, e = F{U) + 0(L _W ), (33) 

where F(U) is a universal function. I 11 Fig. [3] we compare 
the results for the two models: they appear to approach 
the same asymptotic curve with increasing the lattice 
size. Scaling corrections are consistent with the expected 



FIG. 3: (Color online) Plot of R 5 vs U. ACP 2 and 0(8) results 
approach the same asymptotic curve. The dotted horizontal 
and vertical lines correspond to the values I?| ~ 0.5237 and 
U* ~ 1.0383 estimated in the 0(8) vector model. 


L _aj behavior with u a: 0.8. They are much smaller for 
ft < Pc than for /3 > / 3 C . 

The dependence of the data on the inverse temperature 
/3 around the crossing point is consistent with the 0(8) 
results. Fits of the data around the crossing point to the 
first few terms of the expansions appearing in Eq. (1281) 
provide an accurate estimate of the critical point, 

Pc = 4.142(1). (34) 

In Fig.Ejwe show x^ v ~ 2 and versus A', using the esti¬ 
mate (l34l) of p c and the 0(8) estimates (1321) of the critical 
exponents. The data approach asymptotic scaling curves. 
Scaling corrections are larger for Rj, but definitely com¬ 
patible with the expected L _w behavior. Note also (not 
shown) that the scaling curves of R% for the ACP 2 and 
0(8) nicely match after a trivial rescaling of the scaling 
variable (p — /3 C )L 1 / 1 ': if we define Y = X for the ACP 2 
and Y = 3.9X for the 0(8) model, all data fall on the 
same curve when plotted versus Y. Finally, we consider 
the specific heat C v at p c . We expect C v ~ a + cL a l v 
with a/v = 2/u—3 ~ —0.67. The MC data are consistent 
with this behavior. 

In conclusion the numerical analysis of the ACP 2 lat¬ 
tice model provides a robust evidence that its contin¬ 
uous transition belongs to the universality class of the 
0(8) vector model, as predicted by the RG arguments of 
Sec. mi 


IV. RG FLOW FOR N > 4 

In this section we present a FT study of the RG flow 
of the LGW theory E]) with N > 4, i.e. the most gen¬ 
eral <f> 4 theory with traceless hermitian N x N matrix 
fields and parity symmetry. The critical behavior at a 
continuous transition is controlled by the FPs of the RG 
flow, which are determined by the common zeroes of the 
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FIG. 4: (Color online) Scaling behavior of the ratios x/lP v 
and R,£ vs (/3 — /3 C )L R 1 ' with /3 C = 4.142 and the critical 
exponents of the 0(8) universality class, cf. Eq. (1321) . In the 
top figure the dashed horizontal segment indicates the 0(8) 
universal value R% = 0.5237, which is clearly approached by 
the data of the ACP 2 lattice model with increasing L. 


(3 functions associated with the quartic parameters. The 
presence of a stable FP controls the universal features 
of the critical behavior in the case of continuous transi¬ 
tion. The absence of a stable FP implies the absence of a 
corresponding universality class, hence a transition char¬ 
acterized by the same symmetry-breaking pattern must 
be of first order. 


A. The MS perturbative scheme 

We compute the /3 functions of the quartic couplings 
in the MS renormalization scheme [27|, which uses the 
dimensional regularization around four dimensions, and 
the modified minimal-subtraction prescription. Thus the 
RG functions are obtained from the divergences for e = 
4 — d —t 0 appearing in the perturbative expansion of the 
correlation functions of the critical massless theory. 

The procedure is straightforward, see, e.g., Ref. [28| . 
The renormalized couplings are defined from the irre¬ 
ducible four-point correlation function, and the MS p 


functions are 

dv 

, (3 v (u,v) = n— , (35) 

no,t>o U ^ u 0 ,v 0 

where /i is the renormalization energy scale of the MS 
scheme. Here, u and v are the renormalized couplings 
corresponding to zto, vq, defined so that u oc uq/ fj, e and 
v oc v$/if at the lowest order. We compute the /3 func¬ 
tions up to five loops. The complete series for N = 4 are 
reported in App. [Bj 


Q t \ dU 

Pu{u, V) = fl — 


1. One-loop analysis close to four dimensions 


Let us first analyze the one-loop /? functions. They 
read 


Pu 

Pv 


—eu + 


N 2 


27V 2 — 3 TV 2 + 3 2 


37V 


-uv + 


27V 2 


(36) 


—tv + 2uv + 


TV 2 — 9 
37V 


(37) 


The exact normalization of the renormalized variables 
can be easily read from these series. 

Since for TV = 2 and TV = 3 the two quartic terms 
are not independent, an appropriate combination of the 
above P functions must reproduce the p functions of the 
0(3) and 0(8) 4> 4 theories. Indeed, using Eq. (fITl) and 
setting g = u + v/2, we obtain 


Pu + TfPv — Po(3){g) — ~ e 9 + -g -g 2 (38) 

for TV = 2, and 

Pu + 7 ^ Pv = Po(8)(g) — ~ e g + g5 2 (39) 


for TV = 3. These exact relations provide a stringent 
check of the five-loop series for the model ([5]) , which must 
reproduce the cor resp onding series of the 0(3) and 0(8) 
vector model [291, l3Q( for TV = 2 and 3. 

For TV > 4 the FPs of the RG flow are given by 
the common zeros of the p functions © and ©. 
Their stability requires that the eigenvalues of the matrix 
f \j = dPgJdgj (where gi t 2 correspond to u,v) have pos¬ 
itive real part. In the standard e-expansion scheme [3l| . 
the FPs, i.e., the common zeroes of the /3-functions, are 
determined perturbatively as expansions in powers of 
e = 4 — d 1 while exponents are obtained by expanding 
the corresponding RG functions computed at the FP in 
powers of e. 

A straightforward analysis of the one-loop /3 functions 
(157>1) and (1571) finds four different FPs. Two of them have 
v = 0 and are always unstable. We have the trivial Gaus¬ 
sian FP at (u = 0, v = 0), which is always unstable with 
respect to both quartic perturbations. There is also an 
O(M) symmetric FP with M = TV 2 — 1 at 


u = e 


6 


v = 0, 


(40) 


TV 2 + 7’ 











which can be shown, nonperturbatively, to be unstable 
with respect to the operator Tr 'I' 4 . Indeed, such operator 
contains a spin-4 perturbation with respect to the O (M) 
group 0 , which is relevant for any M > 4 to 0(e), and 
for any M > 3 in three dimensions M. 

There are also two FPs with v < 0. One of them is 
stable, the other is unstable. However, they only exist 
for N < N C fi, with 

N cfi = 4=^1 + ^ » 3.506. (41) 

v 2 

For N = N c 0 these two fixed points merge and then, for 
N > N C £, they become complex. For N = 3 the stable 
fixed point merges with the 0(8) fixed point. These re¬ 
sults show that, for integer values of N satisfying N > 4, 
there is no stable fixed point close to four dimensions, 
hence only first-order transitions are allowed. 


2. Five-loop t expansion analysis 

In order to establish the behavior of the system for 
e = 1, we must determine the fate of the stable fixed 
point that exists for N < N c $ close to four dimensions. 
For finite e, we expect a stable and an unstable fixed point 
with v ^ 0 up to N = N c (e). The two fixed points merge 
for N = N c (e) and become complex for N > N c (e). In 
order to compute N c (e ). we expand 

N c (e) = N c , 0 + Y, Nc,ne n , (42) 

n= 1 

and require 

p u (u,v,N c ) = /3 v (u,v,N c ) = 0, 

det Q(u, v , N c ) = 0, (43) 

the last equation being a consequence of the coalescence 
of the two fixed point at N = N c . A straightforward 
calculation gives finally 

N c (e) = 3.5063 - 0.0309e + 0.3229e 2 - 1.2927e 3 

+7.6855e 4 + 0(e 5 ). (44) 

The expansion alternates in sign, as expected for a Borel- 
summable series. Resummations using the Pade-Borel 
method appear to be stable. We obtain N c (e = 1) = 
3.54(1) using the series to order e 3 and N c (e = 1) = 
3.59(2) at order e 4 (the number in parentheses indicates 
how the estimate changes by varying the resummation 
parameters). Apparently, N c varies only slightly as e 
changes from 0 to 1. In particular, this analysis predicts 
the absence of stable FPs for any integer N > 4 in three 
dimensions. 


3. High-order analysis in three dimensions 

Methods based on the expansion around four dimen¬ 
sions allow us to find only the 3D FPs which can be 


defined, by analytic continuation, close to four dimen¬ 
sions. Other FPs, which do not have a 4D counter¬ 
part, cannot be detected. However, the extension of 
this result to the relevant d = 3 dimension fails in some 
cases. For example, this also occurs for the Ginzburg- 
Landau model of superconductors, in which a complex 
scalar field couples to a gauge field: although e-expansion 
calculations do not find a stable FP 33], thus predict¬ 


ing first-order transitions, numerical analyses of 3D sys¬ 
tems described by the Ginzburg-Landau model show that 
they can also undergo continuous transitions, see, e.g., 
Ref. [33, [35| . This implies the presence of a stable FP 
in the 3D Ginzburg-Landau theory- in agreement with 
experiments [36]. Other examples are provided by the 
LGW ‘F 4 theories describing frustrated spin models with 
noncollinear order [13,HI}, the 3 He superfluid transition 
from the normal to the planar phase [39j . and the chi¬ 
ral transitions of the strong interactions in the case the 
U(1) J 4 anomaly effects are suppressed jiol l4lj|. 

Therefore, a more conclusive analysis requires a direct 
study of the 3D flow. This is achieved by an alternative 
analysis of the MS series: the 3D MS scheme without e 
expansion [13, 1H, [Hj . The RG functions /3 U) „ are the MS 
functions. However, e = 4 — d is no longer considered as 
a small quantity, but it is set equal to its physical value 
(e = 1 in our case) before computing the FPs. This 
provides a well defined 3D perturbative scheme which 
allows us to compute universal quantities, without the 
need of expanding around d = 4 [43, |43| . 

We look for stable FPs of the RG flow, with a finite 
attraction domain in the space of the renormalized cou¬ 
plings u and v. The RG trajectories are determined by 
solving the differential equations 

-A^ = /3„KA),u(A)), 
dv 

-A— = /3 v (u( A),u(A)), 
where A € ]0,oo), with the initial conditions 


(45) 


u(0) = v(0) = 0, 
du _ uq 

A —0 ~ M’ 


dv 

~d\ 


= ± 1 , 


A=0 


(46) 


where s parametrizes the different RG trajectories in 
terms of the bare quartic parameters, and the ± sign 
corresponds to the RG flows for positive and negative 
values of vo- In our study of the RG flow we only con¬ 
sider values of the bare couplings which satisfy Eqs. 0 
and ©■ 

The physically relevant results are obtained by resum¬ 
ming the perturbative expansions, which are divergent 
but Borel sunnnable in a large region of the renormalized 
parameters. The resummation can be done exploiting 
methods that take into account their large-order behav¬ 
ior, which is computed by semiclassical (hence, intrinsi¬ 
cally nonperturbative) instanton calculations [251.1281 1321 . 
Relevant results for the large-order behavior of the series 
of the model © are reported in App. [C] 
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B. The 3D MZM perturbative scheme 


In the massive zero-momentum (MZM) scheme mm 
[HI, S3] one performs the perturbative expansion directly 
in three dimensions, in the critical region of the disor¬ 
dered phase, in powers of the zero-momentum renormal¬ 
ized quartic couplings. The theory is renormalized by 
introducing a set of zero-momentum conditions for the 
one-particle irreducible two-point and four-point correla¬ 
tion functions of the 2x2 matrix field 4': 


d2) 


2,bib 2 (P^ ~ (^a 1 b 2 fia 2 b 1 ^ ^a x a 2 £>bib 2 ^ X (47) 


X -Zb, 1 [m 2 +p 2 + 0(p 4 )\ , 


t(4) 

aia2,bib2 ,ciC2,did2 


(0) = Z- 2 m 4 ~ d 


(48) 


X (uU ai a 2 ,bib 2 ,c-LC 2 ,did 2 T ^^a\a 2 ,b\b 2 ,c\c 2 ,d.\d 2 ) ■ 


where U, V are appropriate form factors defined so that 
u oc uo/m and v oc Vo/m at the leading tree order. The 
FPs of the theory are given by the common zeroes of the 
Callan-Symanzik /^-functions 


/3 u {u,v) 


du 
1 dm 


u o,v 0 


Pv{u,v) 


dv 

dm 


■ (49) 

u 0 ,v 0 


The normalization of the zero-momentum quartic vari¬ 
ables u, v is such that their one-loop (3 functions read 


Pu = ~ 


9 47V 4 + 227V 2 - 42 

u + u 2 + 


7V(7V 2 + 7) 2 


■uv 


(50) 


3IV 4 + 307V 2 + 63 


= —v + 


TV 2 (TV 2 

12 


N 2 + 7 


-7) 2 


uv 


27V 4 - AN 2 - 126 

V/ V2 I -\2 


V 2 . (51) 


We compute the MZM perturbative expansions of the 
(3 functions and of the critical exponents up to six loops, 
requiring the computation of 1428 Feynman diagrams. 
The complete expansion for TV = 4 can be found in 
App. [Bj The large-order behaviors of the series are re¬ 
ported in App. O The RG trajectories are obtained by 
solving differential equations analogous to Eqs. (1351) and 
(l46l) . after resumming the /3 functions as outlined in 
App.0 


C. Results 

Some RG trajectories in the renormalized coupling 
space of the LGW theory d5j) for TV = 4 are shown in 
Figs. [5] and [6] for the MS and MZM schemes, respec¬ 
tively, for several values of the ratio s = uo/\vq\. In both 
renormalization schemes most of the RG trajectories flow 
towards the region in which the series are no longer Borel 
summable. In the MZM scheme, for vq < 0 some trajec¬ 
tories flow instead towards infinity. In all cases, we do 



u 

FIG. 5: (Color online) RG flow of the LGW theory © for 
TV = 4, in the MS scheme without e expansion, for several 
values of the ratio s = uo/\vq\ of the bare quartic parameters. 
The curves are obtained by solving Eqs. <Pl5l) with the initial 
conditions (1461) : in the legend we report the value of s and 
the sign of Vo (”+” and correspond to Vo > 0 and Vo < 0, 
respectively). The two solid lines represent the boundary of 
the Borel-summability region, defined by « + »/4 > 0 and 
u + b±v = u + 7v/l2 > 0. 


not have evidence of a stable FP. Analogous results are 
obtained for TV = 6. 

These results imply that there is no universality class 
characterized by the symmetry breakings m and (USD- 
This would suggest a first-order transition. It is also 
possible that more than one transition is present, each 
of them associated to a partial decoupling of some de¬ 
grees of freedom, hence to a different symmetry-breaking 
pattern, as it happens in two-dimensional frustrated XY 
models [45j|. In this case, continuous transition would 
still be possible. 


V. CONCLUSIONS 

We have investigated the nature of the phase transi¬ 
tions in 3D ACP^ -1 models, such as the lattice model 
(CD with J > 0, which are characterized by a global U(7V) 
symmetry and a local U(l) gauge symmetry. 

In order to analyze their critical behavior, we con¬ 
struct the corresponding LGW theory, assuming a stag¬ 
gered gauge-invariant order parameter. This leads to the 
LGW 'F 4 theory ©, where T is a traceless hermitian 
TV x TV matrix, which is symmetric under the global U(7V) 
transformations 4/ —> UA>U\ and the Z 2 transformations 
\F —> — T. For TV = 3, the LGW mode is equivalent to the 
LGW associated with an 8-component real vector field. 
Hence, we predict that, if the ACP 2 model undergoes a 
continuous transition, it should belong to the 0(8) vec¬ 
tor universality class. Note that, at the critical point, 
there is an effective symmetry enlargement, U(3)—»0(8), 
and the same should occur in the low-temperature phase 
as the critical point is approached. The low-temperature 
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u 


FIG. 6: (Color online) RG flow of the LGW theory © for 
N = 4, in the MZM scheme, for several values of the ratio 
s = Mo/| Vo | of the bare quartic parameters. The curves are ob¬ 
tained by solving Eqs. (14511 with the initial conditions (14611 : in 
the legend we report the value of s and the sign of Mo (”+” and 
correspond to Mo > 0 and Mo < 0, respectively). The two 
solid lines represent the boundary of the Borel-summability 
region, defined by u + m/4 > 0 and u + 5 4 m — u + 7m/12 > 0. 


(a) 


(b) 




FIG. 7: (Color online) We draw some lattice configurations 
required by the discussion of Sec. El 


limit. This issue needs further investigation. 


Appendix A: Ground state for the antiferromagnetic 
model with N = 3 


symmetry U(l) x U(2) should be promoted to 0(7). We 
confirm the RG predictions by comparing the FSS be¬ 
havior of the 0(8) vector and ACP 2 models, obtained 
by MC simulations of both lattice models. We note that 
the critical behavior, characterized by the 0(8) critical 
exponents v = 0.85(2) and rj = 0.0276(5), definitely dif¬ 
fers from that of ferromagnetic CP 2 models, for which 
recent studies @, H, Oil |Tl| have provided numerical ev¬ 
idence of continuous transitions with critical exponents 
v = 0.536(13) and rj = 0.23(2). 

In the case of ACP^ -1 lattice models with a higher 
number of components, i.e. N > 4, the identification of 
the order parameter is more complex. If the order pa¬ 
rameter is a staggered gauge-invariant hermitian matrix 
as for N = 2 and N = 3, the associated LGW theory 
is that given in Eq. ©■ To determine the possible exis¬ 
tence of continuous transitions, we study the RG flow in 
perturbation theory. We compute FT perturbative series 
in two different renormalization schemes up to five and 
six loops, respectively. The analysis of the RG flow does 
not provide evidence of stable FPs. This implies that the 
dynamics of the staggered gauge-invariant modes associ¬ 
ated with the NxN hermitian matrices defined in Eq. © 
does not give rise to continuous transitions. In particular, 
we do not expect continuous transitions characterized by 
the symmetry breaking U (N) -A U(l) x U(IV — 1). Thus, 
if the ACP^' 1 lattice model presents transitions with 
this symmetry breaking, they must be first order. Note 
that it is still possible to have continuous transitions if 
they are associated with a different symmetry breaking, 
as it may arise from a partial decoupling of some degrees 
of freedom. Another possible scenario may arise from the 
relevance of further gauge degrees of freedom which are 
not taken into account by the LGW theory, analogously 
to the case of ferromagnetic CP^' 1 model in the large-IV 


In this section we wish to characterize the structure of 
the ground state of the antiferromagnetic model for N = 
3, showing the emergence of a ferromagnetic order on a 
staggered lattice. Let us first consider a lattice plaquette, 
see Fig. ©a), and let us determine the configurations 
of the four spins z 4 , , 23 , and z 4 that minimize the 

energy. Given the global invariance of the model, it is 
not restrictive to assume that Z\ = (1,0,0). Since the 
model is antiferromagnetic, the energy is minimized if 
neighbouring spins are orthogonal, i.e., if Zi ■ Zj = 0 for 
nearest neughbor sites i and j. Therefore, we have 

z 2 = (0,t? 2 ), z 4 = (0,v 4 ), (Al) 

where i >2 and u 4 are two-dimensional complex unit vec¬ 
tors. If we set Z 3 = ( 231 , 1 ) 3 ) and consider the links con¬ 
necting site 4 with its neighbors, we obtain the conditions 

V 2 ■ v 3 = V 4 ■ v 3 = 0. (A2) 

If t >2 ^ e^Vi (</> is an arbitrary phase), this condition 
implies 1)3 = 0. Therefore, discarding an irrelevant phase, 
we obtain z 4 = Z3. If instead If V 2 = e*^i; 4 , eliminating 
an irrelevant phase we obtain z -2 = z 4 . This analysis 
shows, therefore, that in the ground-state configuration 
two opposite spins in the plaquette are identical. Note 
that the result holds only for N = 3. For N > 4, Eq. (IA2I) 
does not imply V 3 = 0 in the generic case. We only 
obtain that Z 3 belongs to (N — 2 )-dimensional subspace 
containing zi. This may leave open the possibility of 
other symmetry breaking patterns for N > 3. 

Let us now assume that we are dealing with a three- 
dimensional system and let us consider a cube, see 
Fig. ©b). We wish now to show that the dominant 
ground-state configurations have the following structure. 
One can have z 4 = Z 3 = Zq = Zg (we properly fix the 
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phases), while z 2 , z 4 , z 5 , Z 7 are arbitrary spins lying 
in the complex two-dimensional plane orthogonal to Z\. 
Of course, the equivalent arrangement with 2 2 = Z 4 = 
z 5 = Z 7 is also possible. Note that this configuration is 
highly degenerate, as four two-dimensional vectors can 
be arbitrarily chosen. 

To prove the previous statement, we consider all other 
possible arrangements that are consistent with the re¬ 
sult obtained for the plaquette and we show that they 
are less degenerate that the one discussed before. Hence 
they are entropically disfavored and become irrelevant in 
the infinite-volume limit. It is not restrictive to assume 
that Zi = z 3 = ( 1 , 0 , 0 ), since two opposite spins of pla¬ 
quette ( 1 , 2 ,3,4) (we report in angular brackets the sites 
belonging to the plaquette) are necessarily identical. Let 
us now consider the plaquette (5, 6 , 7, 8 ). There are three 
different possibilities consistent with the result valid for 
a single plaquette: (i) z 6 = z 8 and z 6 = z 4 ; (ii) z 6 = z 8 
with z 6 ^ e^zi] (iii) z 5 = Z7. We wish now to exclude 
cases (ii) and (iii). In case (ii) the orthogonality condition 
applied to links ( 1 , 2 ) and ( 2 , 6 ) implies 

zi • z 2 = z e ■ z 2 = 0 (A3) 

Since z 6 / e^zi, these two conditions completely specify 
vector z 2 (the phase is of course irrelevant). Repeating 
the argument to the spins at sites 2 ,4,5,7 we end up with 

z 2 = z 4 = z 5 = z 7 (A4) 


It is therefore a configuration of type (i), but now the 
spins are ordered on the complementary sites 2,4,5,7. Let 
us now consider case (iii). Generically, z 2 ^ e*^Z 5 and 
24 7 ^ e^Zrj- This implies that z 6 and z 8 are uniquely 
defined. Therefore, the configuration is defined by spec¬ 
ifying z 2 , z 4 , and Z 5 , in the two-dimensional complex 
space orthogonal to Zi. This configuration is less degen¬ 
erate than that defined at point (i), hence it is irrelevant 
in the infinite-volume limit (it is entropically suppressed). 
If z 2 = e*^Z 5 and z 4 7 ^ e^z 5 , we obtain the same type of 
degeneracy, since now z 4 , Z 5 , and Zq can be chosen. The 
other two cases give the same result. 


The result for a cube extends trivially to the whole 
lattice, proving that in the ground state we observe two 
different symmetry breakings. First, lattice translational 
invariance is broken with the emergence of a staggered 
symmetry. Second, on one of the two sublattices the 
system orders ferromagnetically, breaking the U (3) sym¬ 
metry down to U{ 2). It is important to note that the 
discussion only applies in three dimensions. In two di¬ 
mensions the results does not hold. Indeed, referring to 
Fig. He), nothing forbids in the two-dimensional case a 
configuration with z 4 = 23 , z 4 = Z5, and Z6 7 ^ Z3. 


I 

Appendix B: High-order field-theoretical perturbative expansions 


In this appendix we report the FT perturbative series of the /3 functions used in our RG analysis of Sec. IIVI We 
only report those for N = 4; the perturbative series for other values of N are available on request. 

The five-loop (3 functions of the MS scheme for N = 4 are 


/3 u (u,v) = 
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The six-loop /3 functions of the MZM scheme for N = 4 are 

/3„(m,m) = -u + u 2 - 0.22544283m 3 + 0.10908673m 4 - 0.06576687m 5 + 0.04692261m 6 - 0.03823225m 7 (B3) 

+0. 63043478mm - 0.20303858m 2 m + 0.15749967m 3 m - 0.11634780m 4 m + 0.10631128m 5 m - 0.10058595m 6 m 
+0.15489130m 2 - 0.08970454uo 2 + 0.11442230mV - 0.10401712uV + 0. 12188821m 4 !; 2 - 0.13382427m 5 m 2 
—0. 01961248m 3 + 0.04186104m; 3 - 0.05255068m 2 m 3 + 0.08208912wV - 0.10859850m 4 m 3 + 0.00498971m 4 
—0.01301539mm 4 + 0.03157946wV - 0.05477669mV - 0.00117098m 5 + 0.00642308mm 5 - 0.01679641mV 
+0. 00054799m 6 - 0.00289150mm 6 - 0.00021554m 7 , 

/3 v (vl, v) = -v + 0.52173913mm - 0.20023805m 2 m + 0.06703734m 3 m - 0.05601881m 4 m + 0.03225406m 5 m (B4) 

—0.03176901m 6 m + 0.15217391m 2 - 0.14632780mm 2 + 0.06972703mV - 0.08566328mV + 0.06017577mV 
-0.07408173mV - 0.02133655m 3 + 0.02700191mm 3 - 0.05410921mV + 0.05068529m 3 m 3 - 0.07920335m 4 m 3 
+0.00462087m 4 - 0.01692328mm 4 + 0.02380079mV - 0.04877347mV - 0.00215837m 5 + 0.00591753mm 5 
-0.01769814m 2 m 5 + 0.00059843m 6 - 0.00351442mm 6 - 0.00029807m 7 . 


Appendix C: Summation of the pertubartive series 

Since perturbative expansions are divergent, resumma¬ 
tion methods must be used to obtain meaningful results. 
Given a generic quantity S{u,v) with perturbative ex¬ 
pansion S(u,v) = CijiAv 3 , we consider 

S(xu,xv ) = E Sfc(u, v)x k , (Cl) 

k 

which must be evaluated at x = 1. The expansion m in 
powers of x is resummed by using the conformal-mapping 
method [2§J that exploits the knowledge of the large- 
order behavior of the coefficients, generally given by 

s fc (M, m) ~ k\ [ -A{u , v)] k k b [1 + Oik- 1 )] . (C2) 

The quantity A(u, v ) is related to the singularity t s of 
the Borel transform B(t) that is nearest to the origin: 
t s = —1/A(u, v). The series is Borel summable for x > 0 
if B{t) does not have singularities on the positive real 
axis, and, in particular, if A(u,v) > 0. The large-order 


behavior can be determined using semiclassical compu¬ 
tations, based on the computations of appropriate in- 
stanton configurations |2^, [28[. For even N, these semi¬ 
classical calculations show that the expansion is Borel 
summable when 

m + 6atm> 0, u+— m > 0, (C3) 

where 6 jv is given in Eq. 0. For odd N we obtain anal¬ 
ogously 

m + bjsrv >0, m + cjvm > 0, (C4) 

where cjv is given in Eq. ©■ Note that the conditions for 
Borel summability on the renormalized couplings corre¬ 
spond to the stability conditions 0 and 0 of the bare 
quartic couplings. In this Borel-summability region we 
have for even N 

A(u, v) = — Max (m + 6atm, m + v/N). (C5) 

For odd N , we should replace u + v/N with m + cjvm. Un¬ 
der the additional assumption that the Borel-transform 
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singularities lie only in the negative axis, the conformal- 
mapping method turns the original expansion into a con¬ 
vergent one in the region (IC3I) . Outside, the expansion 
is not Borel summable. 

Analogously one can derive the large-order behavior of 
the MZM scheme, which is again given by Eq. (IC2I) but 
with 

1 32997 

A(u,v) = —^—— Max (u + b^v, u + v/N) (C6) 
N z + 7 


for even N. For odd values of N, u + v/N should be 
replaced with u + c^v. 

We use the conformal mapping method to resum the 
series taking into account what we know about their 
large-order behavior. The method we use is described in 
Refs. |28|,[32|. Resummations depend on two parameters, 
which are optimized in the procedure. The approximants 
we use depend on two parameters a and 6; we use here 
the notations of Refs. p8l. [32|. 
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